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Abstract 

The  yolk,  developed  in  [16,22],  is  a  key  solution  concept  in  the 
Euclidean  spatial  model  as  the  region  of  policies  where  a  dynamic 
voting  game  will  tend  to  reside.  However,  determining  the  yolk  is 
NP-hard  for  arbitrary  dimension.  This  paper  derives  an  algorithm  to 
compute  the  yolk  in  polynomial  time  for  any  fixed  dimension. 


1      introduction 

In  the  Euclidean  spatial  model  [6,3,5,9,34,  e.g.],  each  voter's  most  preferred 
policy  is  represented  by  an  ideal  point  located  in  9?m:  voters  prefer  policies 
(points)  closer  to  their  ideal  points  under  the  Euclidean  norm.  This  is  a 
widely-studied  voting  model,  with  many  applications  (e.g.  [18,31,30,17,36, 
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27,28,29,1]).  In  this  model,  a  classical  equilibrium  or  core  will  generally 
not  exist  ([32,  e.g.]),  and  in  fact  the  space  collapses  into  chaotic  cycles 
[21,33].  Consequently,  much  effort  has  been  made  in  finding  a  satisfactory 
alternative  to  solve  the  equilibrium  problem  in  the  spatial  model.  The 
yolk,  established  in  [16,22],  has  emerged  as  an  important  solution  concept. 
Defined  as  the  smallest  ball  intersecting  all  median  hyperplanes,  it  is  the 
region  of  policies  where  a  voting  game  will  tend  to  stabilize.  The  yolk  is 
also  important  by  virtue  of  its  close  relationships  to  other  solution  and 
evaluation  concepts,  such  as  the  uncovered  set  [24,25,12,22],  the  Pareto  set 
[14],  the  win  set  [12],  and  Shapley-Owen  power  scores  [13]. 

However,  even  determining  whether  the  yolk  radius  is  0  is  co-NP-complete 
in  arbitrary  dimension  [19,2].  Johnson  and  Preparata  [19]  give  a  polyno- 
mial algorithm  which  can  be  used  to  determine  whether  the  yolk  radius  is 
0,  for  any  fixed  dimension.  This  suggests  the  possibility  of  a  polynomial 
algorithm  which  computes  the  yolk,  for  any  fixed  dimension.  Indeed  it 
had  been  thought  ([22,20,14])  that  a  linear  program  in  m  +  1  variables  and 
0^n"7  constraints  would  compute  the  yolk,  but  unfortunately  it  can  fail 
[39]. 

In  this  paper  we  derive  a  polynomial  algorithm  for  arbitrary  fixed  di- 
mension. As  is  typical  with  pseudopolynomial  time  algorithms,  the  time 
complexity  makes  the  algorithm  impractical  for  moderate  or  large  dimen- 
sion. At  present  this  does  not  make  for  a  serious  limitation,  because  nearly 
all  the  empirical  studies  to  date  have  been  in  two  dimensions.  Poole  and 
Rosenthal,  analyzing  19th  and  20th  century  United  States  Congressional 
roll  call  data,  find  that  2  dimensions  have  nearly  the  same  predictive  power 
as  higher  dimensions  [28].  So  in  at  least  some  applications  an  efficient 
algorithm  for  the  two-dimensional  case  may  be  all  that  is  ever  required. 

2      Preliminaries 

Let  us  fix  some  notation.  The  voter  ideal  points  are  a  set  V  C  3ftm,  where 
|V|  =  n. 

For  any  hyperplane  h  in  9?m,  we  denote  the  two  closed  halfspaces  defined 
by  h  as  h+  and  h~ . 

A  hyperplane  h  is  median  (with  respect  to  V)  iff  each  closed  halfspace 


it  defines  contains  at  least  half  the  voters,  that  is  iff 

\h+  nV\>  n/2  and  \h~  nV\>  n/2. 

The  yolk  of  V  is  the  smallest  ball  intersecting  all  median  hyperplanes  (see 
Figure  1). 

A  median  split  of  V  is  defined  to  be  any  pair  of  sets  (5,  T)  such  that: 

•  SuT  =  V; 

•  \S\  >  n/2; 

•  |T|  >  n/2. 

A  hyperplane  {x  :  p  •  x  =  p0}  is  consistent  with  a  median  split  {S,T)  iff 
S  C  h+  C\V  and  T  C  h~  D  V.  The  family  of  all  median  splits  is  denoted  M. 

The  following  Proposition  is  immediate  but  useful: 
Proposition   1:    A  hyperplane  is  median  iff  it  is  consistent  with  some 
median  split  (S,T)  £  M. 

Originally  [16,22]  the  yolk  was  defined  only  for  n  odd.  Koehler[20] 
proposed  extending  it  to  n  even.  In  this  paper  n  may  be  odd  or  even;  the 
algorithm  is  slightly  more  complicated  in  the  latter  case. 

3      what  determines  the  yolk  radius? 

The  main  obstacle  to  an  efficient  algorithm  is  that  there  are  infinitely  many 
median  hyperplanes.  We  would  like  to  pass  from  the  infinite  to  the  finite  by 
extracting  a  small  subset  of  these  which  determine  the  yolk.  That  is,  our 
goal  is  to  identify  a  crucial  "determining"  subset  of  the  median  hyperplanes, 
such  that  the  yolk  can  be  determined  just  from  that  subset.  This  goal  will 
be  realized  by  Theorem  1  in  this  section. 

The  natural  subset  to  try  to  use  for  Theorem  1  would  be  the  set  of 
extremal  median  hyperplanes.  A  hyperplane  is  extremal  iff  it  contains  m 
points  of  V .  (In  the  degenerate  case  we  would  require  the  m  points  to  be 
affinely  independent).  In  m  dimensions,  the  extremal  median  hyperplanes 
are  precisely  the  crucial  subset  to  determine  whether  or  not  the  yolk  radius 
is  0  [19,35].  For  the  case  m  =  2  extremal  median  hyperplanes  are  known  as 


limiting  median  lines  [22,20],  and  indeed  preliminary  algorithmic  efforts  for 
the  two-dimensional  case  ([22,20])  were  based  on  the  plausible  assumption 
that  the  limiting  median  lines  determine  the  yolk. 

Unfortunately,  this  assumption  can  fail,  even  in  two  dimensions  ([39]). 
Thus  we  have  to  augment  the  extremal  set  to  attain  a  determining  set. 
The  set  we  arrive  at  in  this  section  is  infinite  in  one  sense,  but  only  a  little 
bigger  than  the  extremal  set  in  another  sense.  We  defer  the  statement  of 
Theorem  1  to  the  end  of  the  section,  because  the  necessary  terminology  is 
developed  in  the  course  of  its  proof. 

To  develop  the  determining  set,  we  recapitulate  part  of  the  proof  of  yolk 
convergence  in  [40].  There  upper  bounds  are  sought  on  the  radius  of  the 
yolk,  and  the  main  obstacle  is  again  the  infinitude  of  median  hyperplanes. 
The  proof  strategy  is  to  formulate  a  mathematical  program  whose  solution 
is  the  yolk  radius,  and  then  take  advantage  of  certain  "nice"  properties  of 
the  mathematical  program. 

Let  r  denote  the  radius  of  the  yolk  of  the  set  V .  (Since  V  is  static,  we 
suppress  it  as  an  argument.)  For  any  point  x  G  9Rm,  let  r(x)  denote  the 
radius  of  the  smallest  ball  centered  at  x  which  intersects  all  the  median 
hyperplanes.  We  call  r(x)  the  x-centered  yolk.  Obviously  r  =  infzr(x). 

Next  we  develop  a  mathematical  program  for  r(x). 

The  key  to  the  formulation  is  a  kind  of  duality  often  termed  "polarity" 
(see  e.g.  [37,26])  to  characterize  the  set  of  median  hyperplanes.  Any  hy- 
perplane  {y  :  p  •  y  —  p0}  is  specified  by  the  (m  +  l)-tuple  p,po-  If  moreover 
||p||2  =  1  then  p0  is  the  distance  from  that  hyperplane  to  the  origin.  For 
any  particular  median  split,  consistency  of  p,p0  with  that  split  is  simply 
expressed  as  a  set  of  linear  inequalities. 

In  particular,  for  any  median  split  (S,  T)  G  M  and  any  point  z  G  9fm  the 
nonlinear  program  below  finds  the  consistent  median  hyperplane  farthest 
from  z. 

max  \p-  z-  po|  (1) 


TO 


£p?  =1  (2) 

p-Vi  >  p0  Vv,-  G  S  (3) 

P-Vi  <  p0  Vv{  G  T  (4) 

The  first  constraint  (2)  normalizes  the  hyperplane  (p,Po)  and  the  other  two 


linear  inequalities  (3,4)  ensure  it  is  consistent  with  (S, T). 

Now  the  radius  of  the  z-centered  yolk  is  the  distance  from  x  to  the  far- 
thest median  hyperplane.  By  proposition  1  this  is  the  farthest  hyperplane 
consistent  with  some  median  split  in  M.  Therefore,  r(x)  would  occur  at 
the  maximum  of  (l)-(4),  taken  over  M.  When  we  take  this  maximum, 
the  absolute  value  signs  can  be  removed  from  p0  in  the  objective  function 
(1),  because  the  maximizing  over  the  median  split  (T,  S)  is  the  same  as 
minimizing  over  the  median  split  (5, T). 

To  further  simplify  notation,  we  take  x  =  0  in  the  remainder  of  the 
development  of  Theorem  1.  There  is  no  loss  of  generality  since  r(x)  depends 
only  on  the  relative  locations  of  the  ideal  points  V  to  x. 

Therefore,  the  following  family  of  nonlinear  mathematical  programs  de- 
termines the  value  of  r(0),  the  yolk  centered  at  the  origin. 

r(0)  =  max  max. po  s.t.  (5) 

IIpII  =  i  (6) 

p-Vi>  p0  Vvt-  e  S  (7) 

P  •  Vi  <  po  Vt\-  6  T  (8) 

This  formulation  (5-8)  does  not  appear  particularly  attractive.  That  the 
outer  maximum  is  taken  over  the  exponentially  large  set  M  is  not  promising 
at  first  glance. 

In  [40]  the  goal  is  to  find  upper  bounds  on  r(0).  To  this  end,  constraint 
(6)  is  relaxed  to  the  set  of  linear  constraints  —1  <  Pj  <  lVj  =  l,...ra. 
When  this  is  done,  a  lucky  thing  happens:  the  outer  maximization  is  taken 
over  an  exponentially  large  class  of  linear  programs,  but  all  of  the  linear 
programs  share  the  identical  set  of  basic  solutions,  which  is  only  polyno- 
mial^ large.  Then  the  inner  and  outer  maximizations  may  be  replaced  by 
a  single  maximum,  taken  over  all  basic  solutions  which  are  feasible  in  at 
least  one  member  of  the  class  of  linear  programs.  For  (p,  po)  to  be  feasible 
in  at  least  one  of  the  linear  programs,  it  must  be  consistent  with  at  least 
one  of  the  median  splits  in  M.  By  Proposition  1,  checking  this  is  equivalent 
to  checking  if  (p,  p0)  is  median,  which  is  easy  to  do. 

For  the  purposes  of  algorithm  design,  we  require  an  exact  formulation 
rather  than  an  upper  bound.   The  subproblem  of  the  exact  formulation  is 


almost  a  linear  program;  just  one  constraint  is  quadratic.    This  suggests 
that  we  apply  the  Karush-Kuhn-Tucker  conditions  to  (5-8),  and  perhaps  a 
similarly  lucky  simplification  will  obtain. 
The  program  for  r(0)  is: 

max      max/(p)  =  p0  s.t.  (9) 

{S,T)eM 

E,m=iP,2  =  i  (10) 

p  •  Vi  >  po  Vut-  e  s  (11) 

p-  Vi  <  po  V«i  e  T  (12) 

Apply  the  KKT  conditions  to  the  inner  maximization,  using  multipliers 
7T0  for  the  quadratic  constraint  (10)  and  7r,  for  the  linear  constraints  (11,12). 
When  Vi  e  S  D  T,  the  two  constraints  combine  to  yield  an  equality.  So 
a  single  unrestricted  (in  sign)  dual  variable  7rt  can  be  used  for  the  two 
constraints  in  this  case.  The  KKT  conditions  yield: 


X>?  =  i;        (is) 

i=i 

V  Vi  >Po  Vrt-  e  S;  (14) 

p  •  Vi i  <  po  Vvt-  G  T.  (15) 

7rt(p  •  Vi  —  po)  =  0,  Vi  inS  U  T";  (16) 

V/  =  (0,0, ...  ,0, 1)  =  27To(p1,p2,  •  •  •  ,pm,0)  -  X>,(t/,-,  -1).  (17) 

t 

The  first  constraint  (13)  says  p  is  normalized,  ||p||2  =  1.  The  second 
and  third  constraints  (14,15)  say  that  the  hyperplane  (p,Po)  is  consistent 
with  the  median  split  (S,T).  (These  three  are  just  primal  feasibility).  The 
fourth  complementary  slackness  condition  (16)  says  that  the  multiplier  or 
dual  variable  ^  is  zero  for  any  point  i;,-  not  on  the  hyperplane  (p,Po)-  We 
call  the  Vi  that  are  on  the  hyperplane  "binding"  since  their  constraints 
are  binding.  The  last  condition  (17)  is  a  vector  equation  and  breaks  into 
2, 7T,-  =  1  and  27r0p  =  X^7r,f,.  Thus  this  condition  says  that  p  is  an  affine 
combination  of  the  binding  v,,  up  to  scaling  by  the  27r0  factor. 

To  understand  the  geometric  meaning  of  the  KKT  conditions  on  this 
problem,  first  consider  the  simple  case  where  there  is  exactly  one  binding 


point,  say  vk.  Then  p  =  vk/2n0  and  ||p||  =  1  implies  p  =  v*/||vfc||  (and 
also  27T0  =  ||vjb||).  Since  vk  is  binding,  p  •  vk  =  p0  =  |K||.  So  p  is  simply 
the  vector  t;*,  scaled  to  be  a  unit  vector,  and  po  is  the  distance  from  the 
origin  to  the  point  vk.  The  hyperplane  p,p0  is  the  hyperplane  normal  to 
the  vector  from  0  to  vk. 

More  generally  (see  Figures  2  and  3),  suppose  there  are  /  >  2  binding 
Vi  represented  by  the  set  B.  Then  p  =  Y,%eB  *ivi  where  J2%eB  t«  =  1-  Thus 
p,  scaled  by  27r0,  is  in  the  affine  hull  of  the  binding  points  B.  The  scaling 
factor  is  whatever  is  necessary  to  keep  ||p||  =  1.  Also,  p-f,  =  po  Vi  £  B.  For 
notational  convenience,  denote  the  /  points  of  B  as  v\,  •  •  • ,  Vj.  The  preceding 
equation  implies  p-Vi  =  p-v2  =  •  •  ■  =  p-vt.  So  p-(vk  —  Vi)  =  OV/c  =  2, . . . ,/. 
Geometrically,  p  is  the  vector  from  0  to  the  unique  point  closest  to  0  on  the 
affine  hull  of  the  binding  v,-. 

Taken  together,  the  geometric  properties  of  p  imply  that  p  is  the  nor- 
malized normal  vector  to  the  (/  —  l)-dimensional  affine  space  of  the  binding 
points,  in  the  /-dimensional  space  defined  by  the  origin  and  the  points  of  B. 
That  the  points  of  B  are  binding  forces  (p,  p0)  to  be  the  (m—  l)-dimensional 
hyperplane  containing  this  affine  hull  and  normal  to  p  of  course.  Finally, 
the  second  and  third  constraints  require  that  this  hyperplane  (p,po)  be 
consistent  with  the  split  5,  T. 

The  preceding  gives  the  meaning  of  the  KKT  conditions  for  a  single 
split  (5,  T).  Thus  the  KKT  conditions  reduce  consideration  from  an  infinite 
set  to  a  finite  but  exponentially  large  (|.M|)  collection  of  polynomial  sized 
[\B\)  sets  of  hyperplanes.  When  we  take  the  maximum  over  all  median 
splits,  we  are  lucky  as  hoped.  The  only  place  where  (5,  T)  appears  is  in  the 
requirement  that  (p,po)  be  consistent  with  (5,  T).  Taken  over  M,  splits,  this 
is  transformed  into  the  requirement  that  (p,po)  be  consistent  with  at  least 
one  median  split  (5,  T)  £  M.  By  Proposition  1,  this  is  equivalent  to  the 
requirement  that  (p,  po)  be  a  median  hyperplane.  So  the  similarity  among 
the  subproblems,  together  with  Proposition  1,  further  reduces  consideration 
to  a  polynomial  number  of  hyperplanes. 

We  define  the  set  of  determining  hyperplanes  accordingly  (see  Figures 
2  and  3):  let  B  C  V  where  \B\  =  I  <  m.  Let  Aff(B)  denote  the  (/  -  1)- 
dimensional  hyperplane  that  is  the  affine  hull  of  B.  Let  p  be  the  unit 
vector  coincident  with  the  shortest  line  segment  connecting  the  origin  and 
Aff(B),  and  let  (p,po)  be  the  hyperplane  containing  Aff(B).  Then  (p,po) 
is  a  determing  hyperplane.    We  may  refer  to  it  as  the  determining  hyper- 


plane  of  V  specified  by  B  at  0.  When  V  and  0  are  clear  from  context  we 
will  simply  refer  to  it  as  the  determining  hyperplane  specified  by  B. 

If  a  determining  hyperplane  is  median,  it  is  called  a  determining  median 
hyperplane.  We  can  now  state  the  main  result  of  this  section. 

Theorem  1.  The  determining  median  hyperplanes  suffice  to  determine  the 
radius  of  the  O-centered  yolk.        □ 

We  remark  that  after  translation  of  V  Theorem  1  applies  to  the  x- 
centered  yolk  for  any  x.  We  thus  extend  the  definition  of  determining 
hyperplanes  to  the  case  of  arbitrary  x.  It  is  important  to  see  that  the 
determining  hyperplanes  depend  on  x.  The  set  of  determining  hyperplanes 
is  polynomially  bounded  for  any  fixed  x,  although  their  union  is  in  general 
infinite. 

Notice  that  the  determining  hyperplanes  specified  by  sets  B  where 
\B\  =  m,  are  precisely  the  class  of  extremal  hyperplanes  discussed  at  the 
beginning  of  this  section.  Theorem  1  shows  that  these  are  needed  to  deter- 
mine the  yolk,  but  they  must  be  augmented  by  the  hyperplanes  specified 
by  smaller  sets.  The  extremal  hyperplanes  are  the  only  ones  that  are  fixed, 
independent  of  x.  The  other  ones  "swivel"  (or  pivot)  to  be  normal  to  x  as 
x  changes. 

Corollary  1:  The  radius  r(x)  of  the  x-centered  yolk  can  be  determined  in 
polynomial  time  0(nm+1),  for  any  fixed  dimension  m. 

Proof:  The  number  of  determining  hyperplanes  at  x  is  the  number  of 
choices  for  B: 

For  each  choice  of  B,  the  hull  Aff(B),  the  determining  hyperplane  specified 
by  B,  and  its  distance  to  x  obviously  can  all  be  computed  in  0(1)  (constant) 
time  (m  is  fixed).  Moreover,  whether  the  determining  hyperplane  is  median 
can  be  checked  in  0(n)  time.  Thus  the  distance  of  the  farthest  determining 
median  hyperplane  to  x  can  be  computed  in  time  0(nm+1).       □ 


4      a  polynomial  time  algorithm  for  the  yolk 
in  fixed  dimension 

In  this  section  we  develop  a  polynomial  algorithm  to  compute  r.  We  need 

a  couple  of  lemmata  first. 

Lemma  1:  r(x)  is  convex  and  continuous  in  x. 

Proof:    Let  "H  denote  the  set  of  all  median  hyperplanes.    (The  set  V.   has 

infinite  cardinality,  in  general.)  For  any  h  £  M  define  d(x,  h)  as  the  distance 

from  x  to  h.  Then 

r(x)  =  supd(x,/i) 
hey 

Now  d(x,h)  is  convex  in  x  for  all  fixed  /iE)/.  Hence  r(x)  is  convex  since 
the  supremum  of  convex  functions  is  convex.  Continuity  also  follows  easily 
because  the  functions  d(x,  h)  are  a  uniformly  continuous  class.        □ 

We  temporarily  make  a  nondegeneracy  assumption,  that  the  points  V 
are  in  totally  general  position.  Precisely,  we  assume  the  following:  all 
subsets  S  C  V,  \S\  <  m  +  1,  have  full  affine  dimension,  i.e.,  the  points  in 
each  such  S  are  affinely  independent.  This  enables  a  simpler  development 
of  the  algorithm  for  the  nondegenerate  case  (Lemma  2  and  Theorem  1). 
It  will  turn  out  that  the  same  algorithm  works  correctly  in  the  degenerate 
case.  This  extension  will  be  treated  later  (Lemma  3  and  Corollary  2). 

Lemma  2:  Suppose  V  is  nondegenerate  and  fewer  than  m  +  1  determining 
median  hyperplanes  are  binding  at  z.  If  |V|  is  odd  then  z  is  not  the  yolk 
center.  In  the  case  |V|  even,  if  moreover  z  is  not  the  midpoint  of  the  short- 
est line  segment  connecting  two  determining  median  affine  hulls  AfF(51), 
AfT(-B2)  where  \BX\  +  \B2\  <  m  +  1,  then  z  is  not  the  yolk  center. 
Proof:  Let  z  satisfy  the  hypotheses  stated  for  the  case  |V|  odd.  We  wish 
to  define  a  direction  to  move  from  z  along  which  r(x)  is  decreasing.  Let 
J  denote  the  set  of  determining  median  hyperplanes  binding  at  z.  By 
assumption,  \J\  <  m.  Since  |V|  is  odd,  there  are  no  parallel  median  hy- 
perplanes. Thus  by  nondegeneracy  the  hyperplanes  of  J  must  possess  at 
least  one  intersection  point,  say  w.  Then  if  we  move  i  from  z  towards  w, 
we  simultaneously  move  closer  towards  all  the  binding  median  hyperplanes. 
Intuitively,  this  should  decrease  r(x),  the  distance  to  the  farthest  median 
hyperplane. 


If  all  determining  hyperplanes  were  extremal,  this  intuition  would  be 
rigorously  justified  with  no  additional  argument.  However,  from  Theorem 
1  we  know  that  we  must  also  consider  determining  hyperplanes  specified 
by  sets  of  points  B  with  \B\  <  m.  These  hyperplanes  are  more  complicated 
than  the  extremal  ones,  because  they  vary  with  x  (Aff(B)  does  not  vary 
with  x,  but  the  hyperplane  normal  to  the  line  segment  from  x  to  Aff(B) 
does).  In  particular,  we  must  resolve  two  problems.  First,  it  must  be 
possible  to  move  x  a  positive  distance  from  z  without  some  set  B,  whose 
determining  hyperplane  is  not  median  at  z,  suddenly  becoming  median. 
For  if  this  happened,  the  new  median  hyperplane  might  be  far  from  x  and 
r(x)  might  increase.  Second,  suppose  a  hyperplane  in  J  were  non-extremal 
(\B\  <  m).  In  this  case,  the  determining  binding  median  hyperplane  spec- 
ified by  B  rotates  as  we  move  x  from  z  towards  xv.  We  must  be  sure  that 
we  actually  are  moving  closer  to  it,  despite  its  rotation. 

To  deal  with  the  first  of  these  problems,  observe  that  for  any  B,  the  set 
of  points  x,  at  which  the  determining  hyperplane  specified  by  B  is  median, 
is  closed.  Therefore  its  complement  is  open.  Since  z  is  in  the  complement, 
there  is  an  open  ball  containing  z  within  which  the  determining  hyperplane 
specified  by  B  is  not  median.  By  Theorem  1  the  number  of  such  B  is 
finite.  The  intersection  of  a  finite  number  of  open  sets  is  finite.  Therefore, 
for  sufficiently  small  e  >  0,  the  first  problem  does  not  occur  at  the  points 
z  +  e(w  —  z). 

The  second  problem  is  easily  resolved.  Although  the  determining  hy- 
perplane specified  by  B  changes  with  x,  recally  that  the  affine  hull  Aff(B) 
does  not.  Even  when  x  changes,  its  distance  to  B's  hyperplane  is  still  just 
the  distance  to  the  fixed  affine  hull  Aff(B).  Since  B's  hyperplane  is  or- 
thogonal to  z  and  w  is  in  the  hyperplane,  a  sufficiently  small  movement 
towards  w  must  bring  x  closer  to  Aff(B).  This  completes  the  case  when 
|V|  is  odd. 

When  |V|  is  even,  parallel  median  hyperplanes  may  exist.  If  such  a  pair 
is  in  J  then  no  direction  of  improvement  may  exist  since  it  is  not  possible 
to  move  closer  to  both  simultaneously.  Therefore  we  must  supplement  the 
proof  to  deal  with  the  case  of  parallel  determining  median  hyperplanes. 
Fortunately,  this  case  can  only  arise  in  a  very  limited  number  of  situations. 

Suppose  then  that  J  contains  two  distinct  parallel  median  hyperplanes, 
generated  by  JBi  and  B2.  First  we  observe  JB1nB2  =  <j>  since  otherwise  their 
two  affine  hulls  would  intersect  amd  the  parallel  hyperplanes  would  not  be 
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distinct.  Second,  observe  that  \BX\  +  \B2\  <  m  +  1,  for  if  \BX\  +  \B2\  > 
m  +  2  then  d'\m(Aft(Bi))  +  dim(Aff(B1))>  m  and  by  the  nondegeneracy 
assumption  the  two  affine  hulls  would  have  nonempty  intersection. 

Third,  we  claim  that  there  is  a  unique  shortest  line  segment  connecting 
the  affine  hulls  of  B\  and  B2.  For  suppose  otherwise:  let  xy  and  Wz  be  two 
shortest  segments.  We  have 

II*  -  y\\  =  Ik  -  z\\  =    mm    \\Ul  -  u2\\  ?  o, 

u,€A//(B,J 

where  x,w  G  Aff(Bi)  and  t/,2  G  Aff(B2).  The  space  spanned  by  x,w,y, 
and  2  is  at  most  3  dimensional  and  can  be  visualized  easily.  Clearly  if  the 
two  segments  xy  and  Wz  are  not  parallel  then  (x  +  w)/2  is  closer  to  Aff(B2) 
than  z,  a  contradiction.  Thus  x,y,w,z  are  coplanar.  Moreover  they  are 
the  vertices  of  a  rectangle,  or  y  would  not  be  the  point  in  Aff{B2)  closest 
to  x. 

Algebraically,  this  means  x  —  w  =  y  —  z  ^  0.  Recall  that  an  affine 
combination  is  a  linear  combination  in  which  the  sum  of  the  coefficients 
equals  1.  Since  x  and  w  are  affine  combinations  of  the  points  of  £?l5  x  —  w 
can  be  expressed  as  a  nontrivial  linear  combination 

IS!  I 

x-iy  =  ^A,rt;     ^A,=0, 
t=i 

where  rt  £  Bx.    Similarly  y  —  z  can  be  expressed  as  a  nontrivial  linear 

combination 

|Bi|+|B2| 

i=|Bi|+l 

where  r,  £  B2.  The  equality  of  these  two  expressions  implies  the  nontrivial 
linear  equality 

|Bi|+|S2| 

t=i 
where  r,  £  BiU  B2. 

Since  |J5i|  +  |B2|  <  m  +  1,  this  equation  contradicts  the  nondegeneracy 
assumption  of  affine  independence.  This  verifies  the  third  observation,  that 
there  is  a  unique  shortest  line  segment  connecting  the  two  affine  hulls. 
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Fourth  (and  finally),  we  observe  that  z  must  be  the  midpoint  of  this 
shortest  line  segment.  This  follows  because  the  line  segments  from  z  to 
each  affine  hull  are  of  minimum  possible  and  equal  length,  and  are  normal 
to  the  containing  median  hyperplanes. 

Thus  we  have  characterized  the  situations  in  which  two  determining 
median  hyperplanes  from  J  may  be  parallel.  This  completes  the  proof  of 
Lemma  2.       □ 

The  pseudopolynomial  time  algorithm  for  computing  the  yolk  is  based 
on  Theorem  1  and  Lemma  2. 

Theorem  2:  For  any  fixed  m,  the  yolk  of  n  nondegenerate  points  may  be 
computed  in  polynomial  time  0(n^m+1>  ). 

Proof:  The  number  of  determining  hyperplanes  is  0(nm)  by  Theorem  1. 
The  number  of  (m  +  l)-tuples  of  these  is  0(nm(m+1)).  For  each  (m  +  l)-tuple, 
compute  the  set  of  points  equidistant  to  the  corresponding  set  of  Aff(B) 
(this  is  0(1)  since  m  is  fixed).  Let  the  set  of  all  such  points  be  denoted  E. 
Thus  \E\  =  0(nm(m+1)). 

By  Lemma  2,  the  yolk  center  is  in  E.  For  any  point  e  €  E,  the  radius 
r(e)  may  be  computed  in  time  0(nm+1)  by  Corollary  1.  The  point  c  £  E 
that  minimizes^- (e)  :  e  €  Ej'is  the  yolk  center  and  r(c)  is  the  yolk  radius. 
The  time  complexity  is  0{\E\nm+1)  =  o(nm(m+1)+m+1)  =  0(n(m+1)2)  as 
claimed. 

In  the  case  \V\  even,  we  must  also  compute,  for  every  pair  B\,B<l  : 
Bi  n  52  =  <t>]  \B\\  +  \Bi\  <  m  +  1,  the  midpoint  of  the  unique  shortest 
connecting  line  segment.  There  are  0(nm+1)  such  pairs  and  the  midpoints 
must  be  included  in  the  set  of  potential  yolk  centers  E.  The  size  of  E 
remains  0(nm(m+1))  and  the  algorithm  complexity  is  unchanged.        □ 

Now  we  show  that  the  algorithm  of  Theorem  2  applies  to  the  degenerate 
case  as  well.  Lemma  3  shows  that  slight  perturbations  of  V  will  result  in 
only  slight  perturbations  to  the  yolk. 

Lemma  3:  Let  V  be  any  configuration  and  let  V  be  a  nondegenerate 
infintessimal  perturbation  of  V .  Then  the  radius  and  center  of  the  yolks 
of  V  and  V  differ  only  infintessimal ly.  When  \V\  is  even,  the  center  of  the 
yolk  of  V  differs  only  infintessimally  from  one  of  the  yolk  centers  of  V . 
Proof:  Suppose  V  and  V  are  two  configurations  that  differ  infintessimally, 
so  that  \V  —  V\°°  <  e  for  suitably  small  e.  Let  rv  (x)  denote  the  radius  of  the 
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x-centered  yolk  of  configuration  V .  It  is  obvious  that  for  any  x,  rv  (x)  and 
rv(x)  can  only  differ  infintessimally,  since  each  median  hyperplane  moves 
parallel  to  itself  only  infintessimally.  Thus  if  we  slightly  perturb  V  to  get  the 
ideal  points  in  general  position,  the  yolk  radius  is  perturbed  only  slightly. 

However,  this  does  not  yet  assure  that  the  center  of  the  yolk  is  only 
slightly  perturbed  as  well.  Let  C  denote  the  set  of  points  at  which  minz  rv  (x) 
is  attained.  That  is,  C  is  the  set  of  possible  yolk  centers.1  (Note  C  is  con- 
vex.) Now  for  arbitrary  6  >  0  let  C"  denote  {z  :  3y  G  C,  \\z  -  y\\  <  6},  the 
^-neighborhood  of  C.  Let  D  denote  the  closed  set  which  is  the  complement 
of  C".  By  Lemma  1,  the  minimum  of  rv  (x)  on  D  is  attained  at  a  point  q  on 
the  boundary  of  D  with  C".  Let  u  =  rv  (q)  —  rv  (c)  >  0.  Now  select  e  >  0  so 
that  ||rv(x)  -  rK ||°°  <  1//3.  Select  any  ce  C.  Then  for  all  de  L>,  we  have 
r*{d)  >  rv{d)  -  I//3  >  rv{q)  -  u/3  >  rv  [q)  -  2v /3  =  rv  (c)  +  I//3  >  r*(c). 
Therefore  the  minimum  (actually  all  the  minima)  of  rv  (x)  is  (are)  attained 
in  C",  the  ^-neighborhood  of  C.  Thus  for  arbitrary  6  >  0  there  exist  suf- 
ficiently small  perturbations  of  V  that  do  not  perturb  the  yolk  center  by 
more  than  6.  This  justifies  taking  V  to  be  a  nondegenerate  configuration. 
D 

Lemma  3  permits  there  to  be  arbitrarily  close  nondegenerate  configu- 
rations whose  centers  differ  substantially.  But  for  any  nondegenerate  con- 
figuration, there  exists  a  neighborhood  of  configurations  whose  centers  are 
close. 

Example:  Let  V  be  the  degenerate  configuration  of  four  points  as  the 
vertices  of  a  rectangle,  {(0,0),  (2,0),  (0.1),  (2, 1)}.  Then  C  =  (Jl,a)  :  0  < 
a  <  l}.  If  V  is  perturbed  infintessimally  to  nondegenerate  V ,  then  the 
center  of  the  yolk  of  V  will  be  infintessimally  close  to  either  (1,0)  or  (1,  l), 
depending  on  how  the  short  sides  of  the  perturbed  rectangle  slant. 

Now  we  extend  the  algorithm  to  the  degenerate  case. 

Corollary  2:  For  any  fixed  m,  the  yolk  of  n  points  may  be  computed  in 
polynomial  time  0(n^m+1^  ). 

Proof:  Let  V  have  yolk  center  c  and  radius  r.  If  V  were  degenerate,  it 
could  in  principle  be  infintessimally  perturbed  to  nondegenerate  V.  The 
algorithm  of  Theorem  2,  applied  to  V ,  would  compute  center  c  and  radius 
f.  Now  if  the  algorithm  of  Theorem  2  were  applied  to  V ,  it  would  compute 


1]f  the  yolk  is  unique  (it  always  is  when  n  is  odd)  then  C  consists  of  a  single  point. 
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"center"  say,  c'  and  "radius"  r'.  Consider  the  computations  performed 
by  the  algorithm:  they  involve  the  distances  between  the  points  V  and 
hyperplanes  determined  by  these  points.  Moreover  all  computations  are 
confined  to  the  compact  region  defined  by  the  convex  hull  of  V.  Therefore, 
c'  and  r'  must  be  arbitrarily  close  to  c  and  f ,  respectively,  for  small  enough 
perturbation.  But  by  Lemma  3,  the  latter  are  arbitrarily  close  to  c  and  r, 
the  true  yolk  center  and  radius,  respectively.  Thus  el  and  r'  must  in  fact 
be  the  true  values.        □ 

The  algorithm  of  Theorem  2  is  very  easy  to  parallelize,  because  it  simply 
takes  the  minimum  of  a  large  number  of  independent  calculations.  We 
formalize  this  observation  in  a  corollary. 

Corollary  3:  For  any  fixed  dimension,  the  yolk  of  V  may  be  computed  in 
logarithmic  time  with  a  polynomial  number  of  processors. 

Proof:  The  set  of  equidistant  points  E  may  be  computed  independently 
in  parallel  for  each  (m  +  l)-tuple  described  in  Theorem  2.  The  minimum 
of  a  polynomially  bounded  set  of  numbers  may  be  computed  in  logarithmic 
time.  Also,  the  splitting  up  of  the  tuples  requires  only  logarithmic  time. 

It  remains  to  show  that  for  each  equidistant  point  e,  the  radius  of  the 
e-centered  yolk,  r(e),  may  be  computed  in  logarithmic  time.  In  the  proof  of 
Corollary  1,  each  choice  of  B  could  be  processed  independently  in  parallel 
(assuming  sufficiently  many  processors).  The  first  part  of  the  processing  is 
the  computation  of  Af  f(B)  and  its  distance  to  e:  this  requires  O(l)  time 
(in  fixed  dimension).  The  second  part  of  the  processing  is  the  check  to  see 
if  the  determining  hyperplane  is  median  or  not.  This  may  be  accomplished 
in  logarithmic  time  with  n  processors:  each  processor  handles  one  point 
and  checks  in  constant  time  which  halfspace  (possibly  both)  the  point  is 
in;  the  sums  are  tallied  up  in  logarithmic  time. 

The  algorithm  therefore  requires  only  logarithmic  time  and  0(n^m+1^  ) 
processors.        □ 

In  practice,  with  current-day  hardware,  far  less  than  this  full  amount  of 
parallelism  would  be  implementable.  It  is  clear  from  the  proof  of  Corollary 
3  that  a  nearly  perfect  speed-up  could  be  obtained  with  any  small  number 
of  processors. 

Even  for  m  =  2  dimensions,  the  complexity  estimate  of  Theorem  2  is  too 
large  to  ensure  practicality.  In  the  next  section  we  make  some  improvements 
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in  the  algorithm  and  in  the  complexity  estimation,  to  achieve  a  practical 
method  for  problems  of  moderate  size  in  two  dimensions. 

5      improvements  in  performance  and  time  bounds 

In  this  section  we  improve  the  actual  efficiency  and  sharpen  the  analysis 
of  the  algorithm  for  the  two-dimensional  case,  to  achieve  0(n45)  provable 
worst-case  time  complexity.  If  a  conjecture  of  Erdos,  Lovasz,  Simmons,  and 
Straus  holds,  the  worst-case  complexity  reduces  to  o{ns+€)  for  all  e  >  0. 
Moreover,  based  on  empirical  observation,  we  would  expect  the  algorithm 
to  require  0(n3  logn)  time.  We  also  suggest  a  variant  which  could  possibly 
perform  in  sub-cubic  time  (for  2  dimensions),  though  rigorous  bounds  of 
this  quality  remain  an  open  problem. 

The  first  improvement  has  to  do  with  computing  the  set  of  points 
equidistant  to  the  affine  hulls  Aff(B).  When  all  three  hulls  are  lines, 
or  all  points,  there  is  only  one  equidistant  point,  defined  by  the  intersec- 
tion of  two  lines.  But  if  some  of  the  hulls  are  points  and  some  are  lines, 
there  will  be  two  equidistant  points,  defined  by  the  intersection  of  a  line 
and  a  parabola.  In  higher  dimensions  there  may  be  many  (O(l)  for  any 
fixed  m)  equidistant  points.  In  the  previous  section,  we  computed  r(x)  for 
each  of  these  points.  But  this  is  not  necessary.  By  the  convexity  of  d(x,  h) 
in  Lemma  1,  the  distances  to  the  hulls  on  any  convex  combination  of  the 
equidistant  points  are  less  than  or  equal  to  the  convex  combination  of  the 
distances  at  the  equidistant  points.  Therefore,  if  these  hulls  are  the  bind- 
ing (farthest)  from  an  equidistant  point  p,  and  another  equidistant  point  q 
is  closer  to  the  hulls,  then  p  cannot  be  the  yolk  center.  This  means  that 
whenever  there  are  multiple  equidistant  points,  we  may  discard  all  but  the 
closest  of  these. 

Now  consider  a  point  q,  which  we  suppose  is  the  closest  equidistant 
point  to  three  (m  +  1  in  general)  affine  hulls,  at  distance  k.  Suppose  also 
that  the  determining  lines  (hyperplanes  in  general)  are  median.  Then  the 
yolk  radius  must  be  at  least  k,  because  the  triangle  (simplex)  formed  by 
3  (m  +  1)  determining  median  lines  (hyperplanes)  has  an  inscribed  circle 
(ball)  of  radius  k. 

Thus  each  3-tuple  (m-tuple)  of  affine  hulls  provides  a  lower  bound  on 
the  yolk  radius,  if  the  determining  lines  are  median.  In  the  previous  section, 
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we  computed  r(q)  for  each  equidistant  point  q,  and  selected  the  smallest  to 
get  the  yolk  radius.  Instead,  we  discard  all  but  the  closest  in  each  set  of 
equidistant  points;  we  further  discard  any  of  these  whose  containing  lines 
are  not  all  medians;  then  we  select  the  point  among  those  remaining  at 
largest  distance  k  from  the  affine  hulls  it  is  equidistant  to. 

The  revised  algorithm  operates  as  follows:  there  are  0(nm(m+1))  (m+l)- 
tuples  of  affine  hulls.  For  each  of  these,  use  O(l)  time  to  find  the  closest 
equidistant  point  and  compute  the  determining  hyperplanes,  and  O(n)  time 
to  check  if  these  hyperplanes  are  all  median.  Selecting  the  largest  resulting 
distance  takes  no  extra  time.  The  resulting  complexity  is  0(nm  +m+1). 

The  computation  taking  the  most  time  in  the  revised  algorithm  is  the 
check  to  see  if  the  determining  hyperplanes  are  all  median.  If  an  affine  hull 
is  extremal  then  it  is  its  own  determining  hyperplane  (it  doesn't  "swivel" ). 
So  the  determining  hyperplane  is  either  median  or  not,  independent  of 
which  (m  +  l)-tuple  is  being  considered.  It  is  more  efficient  to  precompute, 
once  for  each  extremal  hyperplane,  whether  it  is  median  or  not.  There 
are  only  0(nm(m+1)_1)  (m  +  l)-tuples  of  affine  hulls  in  which  not  all  the 
affine  hulls  are  m-dimensional  (i.e.  extremal,  limiting).  For  these  tuples 
the  processing  time  is  still  0(n);  for  the  0(nm(m+1))-tuples  of  extremal 
hyperplanes  the  precomputations  reduce  the  processing  time  to  O(l).  The 
resulting  complexity  is  0(nm  +m). 

In  the  two-dimensional  case  an  alternative  to  preprocessing  the  extremal 
hyperplanes  is  to  preprocess  the  ideal  points  V  so  as  to  be  able  to  query 
whether  a  line  is  median  in  time  better  than  0(n).  If  0(n2  log  n)  time  and 
0[n2)  space  are  used,  each  such  query  can  be  answered  in  O(logn)  time 
[7].  This  would  result  in  a  complexity  of  0(nm  +mlogn). 

If  the  two  kinds  of  preprocessing  are  both  used,  the  resulting  complexity 
is  still  0(nm2+m),  the  dominant  term  due  to  processing  0(nm2+m)  tuples  of 
extremal  hyperplanes  each  in  constant  time.  So  the  best  2-D  algorithm  we 
have  so  far  has  complexity  0(n6). 

To  further  reduce  the  time  complexity,  notice  that  not  all  of  the  ex- 
tremal hyperplanes  will  be  median.  Computational  experience  makes  one 
suspect  that  there  cannot  be  a  great  many  median  hyperplanes.  In  two 
dimensions  I  have  never  encountered  a  configuration  from  actual  data  with 
more  than  than  2n  distinct  extremal  median  hyperplanes2.  Let  e(V)  denote 

This  is  consistent  with  the  experience  reported  to  me  by  other  researchers. 
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the  number  of  distinct  extremal  median  hyperplanes  of  the  configuration 
V ,  and  let  en  denote  the  maximum  of  e(V)  over  all  nondegenerate  configu- 
rations V  :  \V\  =  n.  The  best  known  upper  bound  on  en  in  two  dimensions 
is  given  by  Erdos,  Lovasz,  Simmons,  and  Strauss: 
Theorem  3  [10]:  en  =  0{n15).       □ 

Erdos  et  al.  [10],  (and  independently  Edelsbrunner  and  Welzl  [8]  3  ) 
construct  configurations  V  with  e(V)  ~  n  log  n.  They  add,  "it  appears 
likely  that  the  lower  bound  obtained  for  en  is  closer  to  the  truth  than  the 
upper  bound  ...  we  conjecture  than  en  =  o(n1+i)  for  all  e  >  0"  (p.  149). 

We  combine  Theorem  3  with  Theorem  2  and  the  algorithm  modifications 
discussed  in  this  section,  to  state  improved  bounds  for  yolk  computations 
in  two  dimensions: 

Corollary  4:  In  two  dimensions  the  yolk  may  be  computed  in  time  0(n4  5). 
If  the  Erdos-Lovasz-Simmons-Strauss  conjecture  holds,  then  the  computa- 
tion may  be  achieved  in  time  o(n3+€)  for  all  e  >  0. 

Proof:  By  Theorem  3  there  are  0(n4B)  3-tuples  of  extremal  median  lines. 
With  the  preprocessing  of  extremal  median  lines  as  discussed,  these  may 
be  processed  in  time  O(l)  each.  There  are  0(n15+1  5+1)  =  0{nA)  3-tuples, 
each  comprised  of  two  extremal  median  lines  and  one  point  (\B\  —  l).  With 
the  preprocessing  of  V  for  halfspace  queries  from  [7],  the  median  queries 
will  require  O(logn)  time  each.  So  these  3-tuples  contribute  0(n4  log  n) 
to  the  total  processing  time.  The  other  kinds  of  3-tuples  will  require  only 
0(n35  log  n)  and  0(n3  log  n)  time  respectively.  Thus  the  total  processing 
time  required  is  0(n45  -f  n4  log  n)  =  0(n45)  as  claimed. 

Now  suppose  the  ELSS  conjecture  is  true.  If  en  =  o(n1+s)  then  there 
are  o(n3+3S)  3-tuples  of  extremal  median  lines  to  be  processed  in  O(l)  time 
each.  There  are  o(n3+2S)  3-tuples  (of  two  lines  and  one  point  each)  to 
be  processed  in  O(logn)  =  o(ns)  time  each.  The  other  kinds  of  3-tuples 
contribute  less.  The  overall  complexity  is  o(n3+e)  where  e  =  36.       □ 

We  remark  that  if  in  practice  e(V)  ~  n,  as  so  far  observed,  then  the 
observed  time  complexity  will  be  0(n3  log  n).  In  this  situation  the  dominant 
term  would  be  from  the  ~  n3  3-tuples  consisting  of  2  extremal  lines  and  one 


3the  upper  bound  of  0{nlb)  in  [8]  on  distinct  feasible  median  splits  does  not  give  an 
upper  bound  on  en  because  many  lines  may  be  consistent  with  the  same  split. 
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point  each.  The  processing  time  would  be  O(logn)  per  3-tuple,  to  check  if 
the  determining  line  of  the  point  is  median. 

For  general  dimension,  finding  good  upper  bounds  on  e(V)  remains  an 
open  extremal  combinatorial  problem.  Erdos  et  al.  remark  (p.  149)  that 
their  upper  bound  does  not  generalize  easily.  Here  we  prove  an  0(nm~  5) 
upper  bound  on  the  expected  value  of  e(V),  in  m  dimensions. 

Proposition  2:  Let  the  n  points  of  configuration  V  be  sampled  indepen- 
dently from  any  nondegenerate  distribution  fi  on  9Jm.  Then  the  expected 
number  of  extremal  median  hyperplanes,  e(V),  is  0(nm~5). 

Proof:  Let  Vi,. . .  ,vn  denote  the  points.  Let  h  be  the  extremal  hyperplane 
specified  by  vi,...,vm.  Suppose  Vi,...,vm  were  sampled  first.  Then  h 
would  be  fixed,  and  we  could  compute  /x(/i+),  the  probability  that  a  point 
sampled  randomly  from  n  falls  in  one  of  the  halfspaces  defined  by  h.  Since 
ix  is  nondegenerate,  fx(h)  =  0.  Now  h  will  turn  out  to  be  median  iff  exactly 
half  the  remaining  n  —  m  points  fall  in  h+ .  The  probability  of  this  event  is 
the  probability  that,  of  n  —  m  identical  independent  Bernoulli  trials,  each 
with  parameter  p  =  n(h+),  exactly  half  are  successful.  This  probability 
is  bounded  by  the  case  p  =  1/2.  In  this  case  the  probability,  from  the 
binomial  distribution,  is  2m~nf "-m)  ~  \j\Jn.  Therefore,  the  probability  h 

is  median  is  at  most  ~  1/y/n. 

By  the  linearity  of  expectation  [11],  the  expected  number  of  median 
extremal  hyperplanes  is  0{nm~5).       □ 

In  the  proof  of  Proposition  2,  if  fi(h)  ^  p,  then  the  probability  h  is 
median  is  smaller  than  the  estimate.  The  slack  in  the  proof  leads  me  to 
conjecture  that  the  bound  of  Proposition  2  holds  in  the  worst  case. 

Applying  Proposition  2  to  the  general  case  gives  0(n^m+1^m~  5))  (m  +  1)- 
tuples  of  extremal  median  hyperplanes.  These  can  be  processed  in  O(l) 
time  each  as  before.  Also  there  are  0(n(m(m~-5)+m~1)  (m  +  l)-tuples  of 
m  extremal  median  hyperplanes  plus  one  determining  hyperplane  specified 
by  B  :  \B\  —  m  —  1.  These  can  be  processed  in  O(logn)  time  each  with 
the  preprocessing  of  V.  The  other  families  of  tuples  are  lower  order  and 
do  not  affect  the  analysis.  The  overall  (expected)  complexity  reduces  to 
0(nm2+m/2--5). 

To  summarize,  the  algorithm  efficiency  can  be  improved  with  appro- 
priate preprocessing.     The  efficiency  also  depends  on  the  magnitude  of 
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e(V).  For  two  dimensions,  the  proven  time  complexity  is  0(n45),  with 
a  conjecture-based  bound  slightly  worse  than  cubic. 

Finally  it  should  be  observed  that  none  of  the  versions  of  the  algorithm 
presented  so  far  has  very  strong  use  of  the  convexity  of  r(x)  from  Lemma  1. 
If  extremal  hyperplanes  sufficed  to  determine  the  yolk,  then  as  McKelvey 
has  observed  [22]  the  yolk  of  V  could  be  calculated  by  solving  an  (m  +  In- 
variable linear  program  in  e(V)  constraints.  The  algorithm  here  performs 
the  analog  of  enumerating  all  or  most  of  the  basic  feasible  solutions.  An 
analogy  with  the  operation  of  the  simplex  method  on  a  min-max  linear 
program  suggests  the  following  algorithm:  (i)  preprocess  extremal  hyper- 
planes and  V;  (ii)  solve  the  linear  program  of  [22]  to  get  a  "hot  start"  initial 
solution  x  equidistant  by  amount  d  to  the  affine  hulls  of  sets  Bx, . . . ,  5m+1; 
(iii)  compute  for  the  initial  solution  the  sum  of  the  violations,  i.e.  the  sum 
over  all  determining  median  hyperplanes  h  at  x,  of  [<i(x,/i)]+;  (iv)  while  the 
sum  of  the  violations  exceeds  0,  iteratively  improve  the  current  solution  by 
exchanging  one  of  the  i?t  for  another  B,  where  improvement  is  measured 
by  the  sum  of  the  violations.  With  0(nm)  choices  for  B,  and  the  compu- 
tation of  the  sum  of  violations  requiring  0(nm),  such  an  algorithm  might 
perform  in  about  0(n2m),  as  compared  with  about  0{nm  )  for  the  enumer- 
ative  algorithm.  In  the  two  dimensional  case  the  algorithm  would  require 
0(n3)  (0(n2+€)  conjectured)  time  per  iteration.  A  theoretical  analysis  of 
this  algorithm,  coupled  with  a  proof  of  the  ELSS  conjecture,  might  at- 
tain a  sub-cubic  time  complexity.  (Notice  a  sublinear  number  of  iterations 
would  be  required.)  While  this  might  not  improve  much  on  n3+£  for  the 
two  dimensional  case,  it  could  make  yolk  computations  in  three  dimensions 
a  practical  possibility. 
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